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1 Introduction 



Quantum electrodynamics is the prototype of a quantum field theory, and its successes in describing 
the interactions of leptons and photons have been spectacular. Nevertheless, continuous quantitative 
tests of QED, particularly at the level of high precision, are important. The positronium system, a 
two-body bound state consisting of an electron and a positron, provides a clean testing ground of QED 
because the effects of the strong and the electroweak interactions are negligible, even at the present 
accuracy of experimental measurements. The existence of positronium was predicted in 1934 j^] based 
on the relativistic quantum theory developed by Dirac and experimentally verified at the beginning of 
the 1950s [||. For the ground state hyperfine splitting, the energy difference between the 1^5*1 (ortho) 
and 1^5*0 (para) states, the most recent experimental values read HI 

W = 203 389.10(74) MHz (1) 

and I, I 

W = 203 387.5(1.6) MHz. (2) 

They represent a precision of 3.6 and 7.9 ppm, respectively, which makes the calculation of all O(a^) 
(NNLO) corrections to the leading and next-to-leading order expression mandatory. Since the dom- 
inant contribution to the hyperfine splitting is of order meO^, NNLO corrections correspond to the 
contributions of order meO®. Including also the known order meO^ln^a"^ contributions |^, 0] the 
theoretical expression for the hyperfine splitting reads|^ 



W = nip a 



12 7rV9 2 y V24 J Sir 



(3) 



where a is the fine-structure constant. At order mga^ it is convenient to distinguish between four 
different sorts of corrections: non-annihilation, single-, two- and three-photon annihilation corrections. 
The two- and three-photon annihilation contributions have been calculated analytically in Refs. 
and [l^], respectively. The single-photon annihilation contributions have recently been determined in 
Refs. [0, In Ref. an analytic result has been presented and in Ref. [0] a numerical one; the 
two results are in agreement. For the non-annihilation contributions, three different results exist in 
the literature [0, 0, |l^, |l^, where Refs. ||l^, 14, |l^ have presented numerical results and Ref. ||l^ 



analytical ones. The results of Refs. [0] and are in agreement. 

A modern and very economical method to calculate non-relativistic bound state problems 
is based on the concept of effective field theories. This approach was first proposed in Ref. ||l^ . 
The effective field theoretical approach to the positronium bound state problem is based on the 
existence of widely separated scales in the positronium system. The physical effects associated to these 
scales are separated by reformulating QED in terms of an effective non-relativistic, non-renormalizable 
Lagrangian, where the low scale effects correspond to an infinite set of operators and the high scale 
effects are encoded in the coefficients of the operators. It is the characteristic feature of the effective 
field theoretical approach that it provides a set of systematic scaling (or power counting) rules that 
allow for an easy identification of all terms that contribute to a certain order in the bound state 



calculation. The results presented in Refs. |12, |l^, 16| have been obtained within an effective field 
theory approach. It is the purpose of this paper to present details of the analytical calculation of 

^ We use natural units, in which Ti — c = 1. The term oc rriea® lna~^ has been determined in Ref. |8|. 
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the order mgO^ single-photon annihilation contribution to the hyperfine splitting presented recently 
in Ref. |2|. 

The program of this paper is as follows: in Sec. |2| we give an overview of the effective field theory 
approach to the positronium bound state problem, and we explain the various steps in the calculation of 
the single-photon annihilation contributions to the hyperfine splitting. Section ^ contains a discussion 
of the subtleties of the cutoff regularization prescription that we use in our calculation. In Sec. ^ we 
describe in detail the calculation of the two-loop short-distance coefficient that is needed to determine 
the single-photon annihilation contributions to the hyperfine splitting at order rriea^, and in Sec. |5|we 
present the bound state calculation, which leads to the final result. A generalization of the result for 
the single-photon annihilation contributions to the hyperfine splitting to general radial excitations is 
given in Sec. ^. Section outlines the status of the theoretical calculations of the hyperfine splitting, 
and Sec. ^ contains a summary. At the end of this work we have attached an appendix where we give 
a collection of integrals that is useful for the matching calculation. 



2 The Conceptual Framework 

The dynamics of a non-relativistic e~^e~ pair bound together in the positronium is governed by three 
widely separated scales: rrig, mgO and mgO'^. Because we are dealing with a Coulombic system, 
where the electron/positron velocity v is of order q ~ q), we could equally well talk about the 
scales nieV and rriev'^ instead of niea and niea'^- These three scales govern different kinds of physical 
processes of the positronium dynamics. The hard scale nie is associated with e^e~ annihilation and 
production processes, the dynamics of the small component and photons with virtuality of the order 
of the electron mass. The soft scale nieV governs the binding of the e~^e~ pair into a bound state 
and directly sets the scale of the size of the bound state wave function, the inverse Bohr radius. The 
ultrasoft scale rrif.v'^ is of the order of the binding energy and governs low virtuality photon radiation 
processes. These processes are associated with higher Fock states, where one has to consider the 
extended system e~^e~'y rather than only an electron-positron pair. Because the interactions between 
the e'^e~ pair associated with a low virtuality photon can arise with a temporal retardation, the effects 
caused by these higher Fock states are called "retardation effects" . The Lamb shift in hydrogen is the 
most famous effect of this sort. The effective field theoretical approach uses the hierarchy of these 
scales (me ^ rriea ^ rrif.a'^) to successively integrate out momenta of the order of the hard and the 
soft scale, and, by the same means, to separate the effects associated with them. In this section we 



give a brief overview onto the conceptual issues involved in this method following Refs. ||l^, 17, 18, |19(| . 
It is the strength of the effective field theoretical approach that it provides systematical momentum 
scaling rules (also called power counting rules) which allow an easy identification of all effects that 
have to be taken into account for a calculation at a specific order. We apply these scaling rules to 
show that retardation effects do not contribute to the hyperfine splitting at order rriea^. 

NRQED is the effective field theory, which is obtained from QED after all hard elec- 
tron/positron and photon momenta, and the respective antiparticle poles associated with the small 
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components have been integrated out. The NRQED Lagrangian reads 



m 



c 
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(4) 



where and x are the electron and positron PauH spinors; Dt and D are the time and space compo- 
nents of the gauge covariant derivative D^, E"^ = F^^ and = ^e^^^F^^ are the electric and magnetic 
components of the photon field strength tensor, and e is the electric charge. The short-distance coef- 
ficients C2,C4,CF,CD,cs,di,d2, which encode the effects from moments of order me, are normalized to 
one at the Born level. The subscripts F, D and S stand for Fermi, Darwin and spin-orbit. In Eq. 
only those terms are displayed explicitly that are relevant to the calculation of the single-photon an- 
nihilation contributions to the hyperfine splitting at order mpU^. The four-fermion operators in the 
last line of Eq. (^) are of particular importance because their coefficients encode the short-distance 
(i.e. hard momentum) effects of the single-photon annihilation process. In general, at order mpa^ for 
the ground state hyperfine splitting, the four-fermion operators shown in Eq. ^ would also contain 
short-distance effects from annihilation into three photons as well as non-annihilation effects^, but 
these are not considered here. The corresponding contributions to the hyperfine splitting have been 
computed elsewhere, using other techniques (see references given in Sees. Q and 0). In the following 
we show explicitly that for the NNLO calculation intended in this work we need the perturbative 
expansion of the constant di at order a^. For ^2 the Born contribution is sufficient. 

From the above Lagrangian, one may derive explicit Feynman rules, after fixing the gauge. As 
is well known, the most efficient gauge for non-relativistic calculations is the Coulomb gauge. In that 
gauge, the Coulomb (or longitudinal) photon (the time component of the vector potential) has an 
energy (i.e. ko) independent propagator, (^o^'^) — This means that the interaction associated 

with the exchange of a Coulomb photon corresponds to an instantaneous potential. The power count- 
ing of diagrams containing instantaneous potentials is particularly simple because an instantaneous 
propagator has no particle pole, i.e. the scale of k is set by the average momentum of the fermions 
~ mv (~ ma in the bound state). On the other hand, the transverse photon (the spatial component 
of the vector potential) has an energy-dependent propagator, of the form (AiA^) ~ {J2 ^i^j)/ i^o ~ 
where the e are the physical (transverse) polarization vectors. In that case, the propagator has a 
particle pole, and and k can be of order mv and also of order mv'^ (with the condition that k^ < \k\ 

^ The effects associated with the two-photon annihilation process would be encoded with the spin singlet operator 
{i'^^2X*) (X^'^2V) (see e.g. Ref. M). 
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([p!9|, 0)). This has the consequence that NRQED diagrams containing transverse photons involve 
contributions from these two scales and, therefore, do not contribute to a unique order in v (or a in 
a bound state) . This fact can be easily illustrated in the context of old-fashioned (or "time-ordered" ) 
perturbation theory in which the integration over the energy components (via the residues) is done 
from the very beginning, and one only has to integrate over the spatial momentum components. Usu- 
ally, the covariant approach is preferred over the old-fashioned perturbation theory because a single 
covariant diagram contains several time-ordered configurations (which are recovered by performing the 
contour integration over the energy components). However, the advantage is lost in a non-relativistic 
application since the different time-ordered diagrams generally scale differently and it is in fact a dis- 
advantage to combine them together. In old-fashioned perturbation theory, one finds that a diagram 
containing an electron-positron pair and a transverse photon will contain a propagator of the form 
(see 1 17] for more details) 

J_ I , (5) 

\k\ Pg^^/me — p^/rue — (p — ky /2me — \k\ ' 

where p^^^t — P — Ta^v are the external and loop momenta of the fermions (we are working in the 
centre-of-mass frame). From this, one can see that different contributions arise depending on whether 
the scale of |fc| is set by p ~ rrieV or by p^ /me — rjieV^. The effects associated with the latter scale 
are the retardation effects. The contributions from both momentum regions will not contribute to the 
same order in v. It is, however, possible to generalize NRQED is such a way that the contributions 
associated to the different scales are coming from separate diagrams. This is achieved by simply 



Taylor-expanding the NRQED diagrams containing Eq. (^) around k ~ nieV and around k ~ nieV 



2 



d^k— = - ^ ' - meV^. (7) 



(the latter expansion is equivalent to a multipole expansion of the vertices) ||T^. One finds that the 
lowest order term of the expansion around k ~ m(,v gives a contribution of order 

T7T T7T - 7 - '^e.v , (6) 

|k| |fe| [rrievy 

whereas the lowest order term of the expansion around k ~ rUeV^ gives 

1 1 imeV^'^'^ 

\^\ Plxt/''^e — P^ /^e + k {nieV^ 

This shows that the dominant contribution from the transverse photon exchange comes from the scale 
k ~ TUeV and that, to leading order, the transverse photon propagator reduces to — (which corre- 
sponds to simply approximating the transverse photon propag ator l/{k'^-k^) by -1/fc^). To leading 
order, the diagrams containing transverse photons are therefore also instantaneous and one recovers 
the simple power counting rules valid for the exchange of a Coulomb photon. At sub-leading order, 
things are more complicated, because both expansions must be taken into account but, fortunately, 
the instantaneous approximation will be sufficient for the present calculation, as will be shown below. 

Because the dominant contribution from the exchange of a transverse photon between an 
electron-positron pair is suppressed by v'^ compared to the dominant contribution from a Coulomb 
photon exchange (see the electron/positron-photon couplings involving the B field in Eq. (^)) all 
interactions at NNLO (i.e. up to order with respect to the Coulomb exchange) can be written as 
a set of simple instantaneous potentials. In momentum space representation they are given by 
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(10) 
(11) 
(12) 



Here, A is a small fictitious photon mass introduced to regularize infrared divergences, Szsp are the 
electron/positron spin operators, and a is the fine structure constant; V^^ is the Breit-Fermi potential 
in the Coulomb gauge, which includes the NNLO relativistic corrections to the Coulomb potential 
from the longitudinal and transverse photon exchange. The potentials V4 and V4dor come from the 
four-fermion operators in Eq. (|^ and account for the single-photon annihilation process at leading 
order and NNLO in the non-relativistic expansion. For convenience we will also count the NNLO 
kinetic energy correction in Eq. ( p!^ ) potential. 

Using the potentials given above, it is straightforward to derive the momentum space equation 
of motion for an off-shell, time-independent e'^e~e^e'^ four-point function in the centre-of-mass frame, 
vahd up to NNLO: 



irip 



Gip,q;s) + J ^Vip,q')G{q,q';s) = {2 nf 6^'\p - q) 



(2vr) 



where 



E = \fs — 2 rrie 

is the centre-of-mass energy relative to the electron-positron threshold and 



V{p,q) 



K;oui(p, q) + Vbt.{p, q) + Vi{p, q) + v^dcXp, q) + sh^^Up^ q) 



(13) 
(14) 

(15) 



The equation of motion (|13| ) is a relativistic extension of the non-relativistic Schrodinger equation of 
the Coulomb problem. Because the potentials V^f, V4, V4dor and lead to ultra-violet divergences, 
it is important to consider Eq. ( |l^ ) in the framework of a consistent regularization scheme. The form 
of the short-distance coefficient di depends on the choice of the regularization scheme. We will come 
back to this issue in Sec. ^. 

One can easily establish simple power counting rules showing that the potentials given above 
are all what is needed for our calculation 0: after factoring out the factors of l/mg that appears 
explicitly in the potentials, the only scale left in diagrams containing the potentials above is the 
inverse Bohr radius (p) ~ meO. In order for the final result to have the dimensions of energy a 
diagram containing any of the potentials shown above will generate one more factor of {p) than there 
are factors of inverse electron mass. If there are n factors of 1/me, the diagram will therefore generate 



a factor {p)"^^^ /ml 



m^a 



n+l 



. This is one source of powers of a. In addition, there are sums over 



We set aside subtleties arising in a cutoff regularization scheme. Those are discussed in Sec. |^. 



6 



intermediate states. Those contain a factor l/{Egxt ~ Eint), which scales hke l/(meQ;^). In order to 
cancel this factor of l/mg, the diagram will also generate a factor of (p), which means that each sum 
over intermediate state brings in another factor of 1/a. Finally, one must multiply by the explicit 
factors of a contained in the NRQED vertices and in the short distance coefficients. As a simple 
illustration of the counting rules, we may consider the Coulomb interaction. The potential contains 
no inverse power of mass (so n = 0) and one explicit factor of a. In first order of perturbation theory 
it therefore contributes to order mga^, which is the same order as the contribution coming from the 
leading order kinetic energy. Adding one more Coulomb potential brings in an extra factor of a from 
the vertices, but this is cancelled by the inverse power of a generated by the sum over intermediate 
states. The Coulomb interaction must therefore be summed up to all orders, as is well known. This 
argument also shows that the Coulomb potential is the only interaction that must be treated exactly, 
as all the other potentials contain at least two powers of inverse mass so that adding one of those 
potentials leads to a contribution of order ma^ (or higher). 

From Eq. (^) we can now derive directly the formula for the single-photon annihilation contri- 
bution to the hyperfine splitting at order niea^, VFnnl'^". Because the para-positronium state does not 
contribute, owing to C invariance, one starts with the well-known n = 1, ^Si e'^e^ bound state wave 
function of the non-relativistic Coulomb problem and determines M^nnlo"" via Rayleigh-Schrodinger 
time-independent perturbation theory. Because we are interested in the single-photon annihilation 
contributions only corrections with at least one insertion of V4 or V4dcr have to be taken into account. 
The formula for W^~^ at order mgO^ then reads 

VF'"^'"" = {l^Si\V4\l^Si) + {l^Si\V4i,,\l^Si) 



+ 



Eq — El 

_l){l 



^ ^0 — -t^i 



+ ..., (16) 



where 1 1 ) represent normalized (bound state and continuum) eigenstates to the Coulomb Schrodinger 
equation with the eigenvalues Ef, \ l^Si ) and Eq = —rrieO? /A denote the state and binding energy of 
the n = 1, ^Si Coulomb bound state. Using the counting rules developed above it is easy to show 
that Eq. ( |T^ ) is all we need to determine the ground state hyperfine splitting to order mgO^: the 
four-fermion operator V4 contains two powers of inverse mass and one explicit factor of a (with the 
Born level value for the coefficient di, see Eq. (p!oD). The contribution of this interaction is therefore 
of order mgO^. In order to obtain the 0{mf.a^) contribution that we are looking for, we therefore 
need to match the coefficient di to two loops, as mentioned above. The operator T^dor, on the other 
hand, contains four powers of inverse mass and therefore contributes already to order nieOi^ with the 



Born level coefficient given in Eq. (11). It is easy to verify that the terms evaluated in second order 
of perturbation theory also contribute to this order if one uses the Born level coefficients in all the 
potentials. Consider for example the term with two insertions of the potential V4. Since there are 
four explicit powers of inverse mass, two explicit factors of a (with di set to 1), and one sum over 
intermediate states, the final contribution is of order mea^"*"^"^ = mgO^. The Breit potential obviously 
contributes to the same order. The operator does not contain any factor of a, but it contains one 
more power of inverse mass and therefore also contributes to order mf>a^ . All other potentials built 
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from the NRQED Feynman rules have higher powers of inverse mass and will therefore be suppressed. 
We note again that the Breit-Fermi potential V^p contains contributions arising from the exchange of 
Coulomb photons and of transverse photons in the instantaneous approximation (i.e. without any ko- 
dependence in the propagator). Since, as we have shown before, the latter contribute already to order 
rUeOe^, we do not need to consider any sub-leading terms coming from the expansions around k ~ rrieV. 
Terms from the expansion around k ~ mev'^ do not need to be considered at all. The instantaneous 
approximation for the transverse photons is therefore sufficient for the present calculation. 

From the above discussion, it is clear that the calculation of W^'^i^q'^ proceeds in two basic 

steps. 

1. Matching calculation - Calculation of the 0{a) and 0{a'^) contributions to the constant di by 
matching the QED amplitudes for the elastic s-channel scattering of free and on-shell electrons 
and positrons via a single photon, e^e~ ^ 7 — > e^e~ , close to threshold up to two loops and to 
NNLO in the velocity of the electrons and positrons in the centre-of-mass frame. This is possible 
because the short-distance effects encoded in di do not depend on the kinematic situation to 
which the NRQED Lagrangian is applied. 

2. Bound state calculation - Calculation of formula on the RHS of Eq. (|l6|). 

The details of the calculations involved in steps 1 and 2 are presented in Sees. ^ and |5[ respectively. 

To conclude this section we would also like to briefly mention a formal way to establish the 
multipole expansion and the counting rules presented above. This is achieved by integrating out 
NRQED electron/positron and photon momenta of order niea. The resulting effective theory has 
been called "potential NRQED" (PNRQED) g. The basic ingredient to construct PNRQED is to 
identify the relevant momentum regions of the electron/positron and photon field in the NRQED 
Lagrangian (^). These momentum regions have been found in Ref. [^]. Because NRQED is not 
Lorentz-covariant, the time and spatial components of the momenta are independent, which means 
that the time and spatial components can have a different scaling behaviour. The relevant momentum 
regions are "soft"^ {k^ ~ nieV, k ~ ruev), "potential" (A;" ~ mgf^, k ~ rngv) and "ultrasoft" {k^ ~ 
nieV^, k ~ ruev'^). It can be shown that electron, positrons and photons can have soft and potential 
momenta, but that only photons can have ultrasoft momenta. A momentum region with fcO 
k ~ niev'^ does not exist. PNRQED is constructed by integrating out "soft" electrons/positrons and 
photons and "potential" photons. In addition, the "potential" photon momenta have to be expanded 
in terms of their time component, because the latter scales with an additional power of a with respect 
to the spatial components. The exchange of "potential" photons between the electron and the positron 
then leads to spatially non-local, but temporally instantaneous, four-fermion operators that represent 
an instantaneous coupling of an electron-positron pair separated by a distance of order the inverse 
Bohr radius ~ rriea. The coefficients of these operators are a generalization of the notion of an 
instantaneous potential. Generically the PNRQED Lagrangian has the form 

/:PNRQED = ^NRQED + j (V' {v) V {v) (x^) (0) + . . . , (17) 

where the tilde above >Cnrqed on the RHS of Eq. ( p!7| ) indicates that the corresponding operators 
only describe potential electrons/positrons and ultrasoft photonic degrees of freedom and that an 

* The soft momentum regime has not been taken into account in the arguments employed in Ref. jl7j. However, this 
does not affect any conclusions concerning the ground state hyperfine splitting at order me a®. 
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expansion in momentum components ~ mf.a'^ is understood. To NNLO, the contributions to V are 
just given in Eqs. ^ to (lTl|). Using the scahng of "potential" electron/positron momenta, we see that 
the Coulomb potential scales like mga^, i.e. it is of the same order as the electron/positron kinetic 
energy. Thus, the Coulomb potential has to be treated exactly rather than perturbatively. Prom the 
PNRQED Lagrangian it is straightforward to derive the momentum space equation of motion of an 
off-shell, time-independent {e~^e~){e~^e~) four-point function in the centre-of-mass frame valid up to 
order a^, Eq. (0). Using the momentum scaling rules of PNRQED one can show that retardation 
effects cannot contribute to W^~'^ at order rriea^. Retardation effects are caused by the ultrasoft 
photons, because their low virtuality propagation can develop a pole for the momenta available in 
the positronium system. Choosing again the Coulomb gauge for our argumentation, where the time 
component of the Coulomb photon vanishes, only the transverse photon needs to be considered as 
ultrasoft|. Thus the emission and subsequent absorption of an ultrasoft photon between the electron- 
positron pair are already suppressed by v'^ ~ with respect to the Coulomb interaction owing to the 
coupling of transverse photon to electrons/positrons. To see that an additional power of a arises from 
the corresponding loop integration over the ultrasoft photon momentum, let us compare the scaling of 
the product of the integration measure and the photon propagator in the potential and the ultrasoft 
momentum regime. In the ultrasoft case the product of the integration measure d^k and the photon 
propagator counts as x = a^, whereas in the potential case the result reads x = a^. 
Thus the exchange of an ultrasoft photon is suppressed by an addition power of a with respect to 
the effects of the Breit-Fermi potential (^. In other words, retardation effects cannot contribute to 
^1-7 ann Qrdcr TTT-ga^. We would like to note that PNRQED is designed as a complete field theory 
capable of describing the dynamics of a bound electron-positron pair and ultrasoft photons. Although 
useful for establishing consistent counting rules, its full strength only develops if one explicitly considers 
the dynamics of ultrasoft photons. For cases where the instantaneous approximation is sufficient - 
such as the ground state hyperfine splitting at order mga^ - the introduction of PNRQED is not 
essential. 



3 The Regularization Scheme 

All equations in the previous section have to be considered within the framework of a consistent UV 
regularization scheme. In general, the form of the short-distance coefficients of the NRQED^ operators 
depends on the choice of the regularization scheme. In this work we use a cutoff prescription to 
regularize the UV divergences, where the cutoff A is considered much larger than mga. The infrared 
divergences, which arise in the intermediate steps of the matching calculation to determine the higher 
order contributions to di, are regularized by a small fictitious photon mass A, see Eqs. and (^). 
The use of a cutoff regularization involves a number of subtleties that shall be briefly discussed in this 
section. 

^ The argument is true in any gauge after gauge cancellations. The argumentation is, however, most transparent in 
the Coulomb gauge. 

^ In what follows, when using the notion "NRQED", we actually mean the generalized NRQED or PNRQED, as 
discussed in Sec. hi 
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Figure 1: Routing convention for loop momenta in ladder diagrams. Half of the external centre-of- 
mass energy is flowing through each of the electron and positron lines. 



It is well known that the use of a cutoff regularization scheme leads to terms that violate gauge 
invariance and Ward identities. These effects, however, are generated at the cutoff and are, therefore, 
cancelled by corresponding terms with a different sign in the short-distance coefficients of the NRQED 
operators. Thus, gauge invariance and Ward identities are restored to the order at which the matching 
calculation has been carried out. Another subtlety is that a cutoff scheme is only well defined after a 
specific momentum routing convention is adopted for loop diagrams in the effective field theory. It is 
natural to choose the routing convention employed in the equation of motion (13). For clarity we have 
illustrated this convention in Fig. |l|. Because we only need to consider interactions in our calculation 
that are instantaneous in time, only ladder-type diagrams have to be taken into account. 

A very important feature of a cutoff regularization scheme is that it inevitably leads to power- 
counting-breaking effects. This means that NRQED operators can lead to effects that are below 
the order indicated by the momentum scaling rules described in the previous section. Examples 
of this feature will be visible in the matching, and the bound state calculations presented in the 
next two sections. These power-counting-breaking effects are a consequence of the fact that a cutoff 
regularization does not suppress divergences of scaleless integrals (as do analytic regularization schemes 
like MS). To illustrate the problem let us consider the third term on the RHS of Eq. (16). This term 
contains the contribution to W^''' coming from two insertions of the four-fermion operator V4. 
According to the momentum scaling rules described in the previous section, it can only contribute at 
order mgO^. However, the bound state diagram with two V4 operators is linearly divergent and, using 
the momentum cutoff A, the result has, for dimensional reasons, the form AnieCt^ + BAa^ where A 
and B are finite constants (modulo logarithmic terms). If one counts the cutoff to be of the order 
of nie, then the diagram would contribute to both orders mgO^ and mgO^. Even worse, by including 
sufficiently high order operators (in the p/nie expansion), one can easily convince oneself that an 
infinite number of operators, having much higher dimension than indicated by the counting rules, 
would contribute to any given order in a, starting at order 171^0^. However, contributions coming 
from those higher-dimension operators can only arise in the form of explicit cutoff-dependent terms 
and not as constants. As for the effects that violate gauge invariance and Ward identities, all terms 
depending on the cutoff are cancelled in the combination of the bound state integrations and the 
short-distance coefficient. In our perturbative calculation we can therefore simply ignore that the 
scaling violating terms coming from operators with dimensions higher than indicated by the counting 
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1 2 

Figure 2: Graphical representation of the NRQED single-photon annihilation scattering diagrams at 
the Born level and at NNLO in the non-relativistic expansion. 

rules exist. For our calculation this means that we only have to calculate the terms presented on the 
RHS of Eq. (|l6|), excluding the higher order corrections of di in the third and fourth terms. 

Finally, we would like to mention that we implement our cutoff regularization scheme in such 
a way that only divergent integrations are actually cut off. This choice simplifies the calculations, 
because we can use the known analytic solutions of the non-relativistic Coulomb problem for the 1^5*1 
wave function and the Green function in our perturbative calculation. The fact that we use a specific 
routing convention ensures that this does not lead to inconsistencies. In addition, we impose the cutoff 
only on the spatial components of the loop momenta. 



4 The Matching Calculation 

The single-photon annihilation contributions of the short-distance coefficient of the operator 
{ijj'^ C7(J2X*) {x^ ^2'^'^) are obtained by matching the amplitude for elastic s-channel scattering of an 
e^e~ pair via a virtual photon in full QED in the kinematical regime close to threshold to the same 
amplitude determined in the non-relativistic effective theory NRQED. To determine the short-distance 
coefficient di to order we have to carry out the matching at the two-loop level, including all effects 
up to NNLO in the non-relativistic expansion. Because we regulate infrared divergences in the effective 
theory using a small fictitious photon mass, we have to do the same in the full QED calculation. 

To obtain the single-photon scattering amplitude in NRQED we have to calculate the Feyn- 
man diagrams depicted in Figs. |2|, ^ and |^, where the various symbols are defined in Fig. The 
Feynman diagrams for the single-photon annihilation contributions of the NRQED scattering ampli- 
tude contain, as the formula for W^'"' """" in Eq. (0), at least one insertion of V4 or V4der- As shown 
in the wave equation (|13|), we can use time- independent electron-positron propagators also for the 
NRQED scattering amplitude. This is possible because all interactions are instantaneous in time, i.e. 
the loop integration over the energy components of the NRQED electron and positron propagators is 
trivial by residue (see Eq. (|2^)). Because the single-photon annihilation process is only possible for the 
electron-positron pair in a ^Si spin triplet state, we only need to consider the Breit-Fermi potential 
in the ^Si configuration: 
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Figure 3: Graphical representation of the NRQED single-photon annihilation scattering diagrams at 
the one-loop level and at NNLO in the non-relativistic expansion. 



where the angular integration it carried out over the angle between p and q. We have eliminated 
the photon mass in the first term on the RHS of Eq. (18) because it does not lead to any infrared 
divergences. An additional simplification for the NRQED calculation is obtained by replacing the 
centre-of-mass energy relative to the e'^e~ threshold, E = ^/s — 2me, by the new energy parameter 
Po, which is defined as 

(19) 



4 



The parameter pQ is equal to the relativistic centre-of-mass three-momentum of the electron/positron 
in the scattering process. At NNLO in the non-relativistic expansion we have the relation 



E 



me 
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Po 



+ ... 



(20) 



where we keep the term p^jme in the LO non-relativistic electron-positron propagator. In this con- 
vention, the NNLO kinetic energy correction reads 
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and simplifies the form of an insertion of the kinetic energy correction 
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We emphasize that the introduction of the parameter po is just a technical trick for the matching 
calculation, which does not affect the form of di. Using the spin average 

\ \{^P^aa2X*){x^^2CTi^)] = 2 (23) 

J=1,0 -1 

for the four-fermion operators V4 and V4dcr in Eq- (§)> we arrive at the following results for the diagrams 
displayed in Figs, g | and | {D{k) = me/{k'^ - Po " ^^))- 
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Figure 4: Graphical representation of the NRQED single-photon annihilation scattering diagrams at 
the two-loop level and at NNLO in the non-relativistic expansion. 
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Figure 5: Symbols describing the NRQED potentials that have to be taken into account for the 
matching calculation at NNLO in the non-relativistic expansion. 
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The upper index of the functions /j corresponds to the power of the fine structure constant of the 
diagrams and the lower index to the numeration given in Figs. ^, ^ and Combinatorial factors are 
taken into account. We note that all the above results have been given in the limit A <C po ^ i^e 
and that only the powers of po relevant at NNLO have been kept. A collection of integrals that have 
been useful in determining the results given above is presented in Appendix. The full NRQED 
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amplitude reads 
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The elastic single-photon annihilation amplitude in full QED can be obtained from the elec- 
tromagnetic form factors, which parametrize the radiative corrections to the electromagnetic vertex, 
and the photon vacuum polarization function. The electromagnetic form factors Fi (Dirac) and F2 
(Pauli) are defined through 
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for the e^e~ production vertex, where q = p -\- p' and 0"^,^ = | [7^,7i/]- We need the form factors 
in the limit A ^ po ^ "^e- The one- loop contributions have been known for a long time for all 
momenta [22, 23 1, whereas the two- loop contributions have been calculated in the desired limit in 
Ref. II^. Parametrizing the loop corrections to the form factors as 

- i + o^^V) + ©Vf)(,^) + ..., 
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and using the energy parameter po (Eq. ([19D) the results for the form factors in the threshold region 
at NNLO in the non-relativistic expansion read (A ^ po ^ w-e)- 
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The photon vacuum polarization function 11 is defined as 



(52) 



where j'^ is the electromagnetic current. The one- and two-loop contributions to 11 are also known from 
Refs. [p^, |23] for all values of q^ and read, expanded up to NNLO in the non-relativistic expansion, 
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Including all effects up to NNLO in pQ/nie and taking the "^Si spin average, the QED amplitude reads 
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The short-distance coefficient di is determined by requiring equality of all the terms up to order 
and NNLO in pQ/m^ in the QED and the NRQED amplitudes in Eqs. ( |56|) and (^3|). The result for 
di reads 
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5 The Bound State Calculation 

In the final step we have to evaluate the RHS of Eq. (p!6|). It is convenient to perform this calculation 
also in momentum space representation. In this representation the normalized n = 1, ^Si positronium 
bound state wave function in the non-relativistic limit reads 

MP)^{P\1'S,) = ^^^, (60) 

where 

7 = (61) 
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We also need a momentum space expression for the sum over intermediate Coulomb states in the 
third and fourth terms on the RHS of Eq. (|T6|). This sum is just the Green function of the non- 
relativistic Coulomb problem, where the n = 1, ^5*1 ground state pole is subtracted and E = Eq = 
—rriea'^ /A. A compact momentum space integral representation for the Coulomb Green function has 
been determined by Schwinger [p5|: 
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Taking the limit E ^ Eq the third term on the RHS of Eq. ( p^ ) develops the n = 1, ^Si pole, 
|0o(p)P/(-E'O ~ E). After subtraction of this pole, one finds that the momentum space representation 
of the sum over intermediate states in Eq. (1T6| ) can be written as: 

\l){l\ . . (2 7r)3 5(3)(p-q)m, 



Ei + ie 



P 



•2 — nip El 



- ^e 



nip 



Att a 



nip 



m En 



qY q^ — niEi 



R{p,q) 



(64) 



where 



R{p,q) 



64 7r7^ 



Q,(p2 +^2)2 (g2 -\-j2yz 



(p2 + -y2) (g2 ^ ^2 



472 



p2 + 72 



' + - In ^ 
2 



g2 _^ ^2 



+ 



2^ - 1 
V4>1- 1 



arctan 



4A- 1 



472 (p 



(65) 



(66) 



Details of the derivation of expression ( p5| ) can be found in Ref. |26]. In Eq. (p^), the three terms 
correspond to no Coulomb, one Coulomb, and two and more Coulomb potentials in the intermediate 
state. In the bound state calculation of Eq.(|l^), the no-Coulomb contributions are linearly divergent, 
the one-Coulomb contributions are logarithmically divergent and the R-term contributions are finite. 
In the case of the bound state contributions involving the R terms, the following relation is quite 
useful {q = \q\ 
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The results for the individual contributions of the RHS of Eq. (16) read 
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where, in the last three terms, the resuhs from the no-Coulomb, one-Coulomb and R-terms have been 
presented in separate brackets. We note that for the bound state calculation we have adopted the 
usual energy definition as given in Eq. ([l3|). A collection of integrals, which were useful to obtain the 
results given above, can be found in Ref. 0. 

Adding all terms together and taking into account the corrections to di shown in Eq. (57), we 
arrive at the final result for the single photon annihilation contributions to the hyperfine splitting at 
NNLO Pi 
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The lna~^ term in Eq. (73) was already known and is included in the lna~^ contribution quoted in 
Eq. (P). The single photon annihilation contribution to the constant K defined in Eq. (^ corresponds 
to a contribution of —2.34 MHz to the theoretical prediction of the hyperfine splitting. The same 
result has been obtained in Ref. |11] using the Bethe-Salpeter formalism and numerical methods. 
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6 Radial Excitations 



From the results presented in the previous sections it would be straightforward to determine the 
single photon annihilation contributions to the hyperfine splitting for arbitrary radial excitations n: 
we "only" have to redo the bound state calculation shown in the previous section using the general 
n^Si wave functions and the momentum space Coulomb Green function, Eq. (^), where the n^Si 
pole is subtracted. The matching calculation (i.e. the form of the short-distance coefficient di), which 
does not depend on the non-relativistic dynamics, would remain unchanged. Although such a strategy 
would be perfectly suited to the general spirit of this work, it would be quite a cumbersome task to 
work our all the formulae for arbitrary integer values of n. (Unlike the contributions to the hyperfine 
splitting from the annihilation into two ^] and three photons which are pure short-distance 
corrections and therefore have a trivial dependence on the value of n (oc |(^„(0)p ~ 1/n^), the single 
photon annihilation contributions have a more complicated dependence on the value of n because they 
involve a non-trivial mixing of bound state and short-distance dynamics.) Thus for the calculation 
of the single photon annihilation contributions to the hyperfine splitting for arbitrary values of n we 
use a much simpler method, which one might almost call a "back of the envelope" calculation ||T^ . 
However, we emphasize that this method relies more on physical intuition and a careful inspection 
of the ingredients needed for this specific calculation than on a systematic approach that could be 
generally used for other problems as well. Nevertheless, this method leads to the correct result, and 
we therefore present it here as well, following Ref. jl2[ |. 

We start from the formula for the single photon annihilation contributions to the hyperfine 



splitting, Eq. (|16D, generalized to arbitrary values of n. (This amounts to replacing "I'^S'i" by "n-^^i" 
everywhere in Eq. (|T^). In the following we simply refer to this generalized equation as "Eq. ([l6|)" .) 
Because we use Eq. ([T^ ) only to identify two physical quantities from which the single photon con- 
tributions to the hyperfine splitting can be derived, we can consider the operators as unrenormalized 



objects. The relevant physical quantities are easily found if one considers Eq. (16) in configuration 
space representation, where the operator V4 corresponds to a (5-function. From this we see that Eq. ([T6|) 
depends entirely on the zero-distance Coulomb Green function An = {0\ Ei^-eI I ^ ) (where the 

n'^Si bound state pole is subtracted) and on the rate for the annihilation of an n'^Si bound state 



into a single photon, P„ = ( n | ^4 | n ) -|- [( n 1 14 eI-Ei (^bf + ^-f^kin) I ) + h.c] + ( n | Vi^^ 

(where the effects from Vbf and SH^^i,^ are included in the form of corrections to the wave function). 
Because we have only considered unrenormalized operators. An and Pn are still UV-divergent from 
the integration over the high energy modes. In the NRQED approach worked out in the previous 
sections the renormalization was achieved at the level of the operators. Now, the renormalization 
will be carried out by relating An and P„ to physical (and finite) quantities, which incorporate the 
proper short-distance physics from the one photon annihilation process. For An this physical quantity 
is just the QED vacuum polarization function in the non-relativistic limit and for P„ the Abelian 
contribution of the NNLO expression for the leptonic decay width of a super-heavy quark-antiquark 
n^Si bound state p^]. Both quantities have been determined in Refs. |2^. From the results given 
in Refs. |28, 29 1 it is straightforward to derive the renormalized versions of An and Pn, 
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where Ci = 2^(-3 + f Ca) - + | In 2 and C2 = ;^(^ - Ca) + | In 2 - f|. Here, 7 is the Euler 
constant and ^' the diganima function. Inserting now A^^"" and P^''^" back into expression we 
arrive at 



which leads to IT 



I-7 ann 



1 



2 Q vr 



2 a vr 



^phys 



mt 



(76) 



aim 



11 rripa 



+ 



4n3 



2^1 + ^2 + ^ 



37 1 

24 n2 6 



1 ^ /2nx 
-+ln — 

71 V a / 



7 - ^{n) 



(77) 



For the ground state n = 1 Eq. ( |77] ) reduces to the result shown in Eq. (|73|). Equation (|77| ) has also 



been confirmed by an independent calculation in Ref. |3C] 



7 Discussion 

In Table Q we have summarized the status of the theoretical calculations to the positronium ground 
state hyperfine splitting, including our own result. To order mga^ the contribution that it logarithmic 
in a and the constant one are given separately. The constant terms are further subdivided into non- 
recoil, recoil, and one-, two- and three-photon annihilation contributions. The non-recoil corrections 
correspond to diagrams in which one or two photons are emitted and absorbed by the same lepton. 
(One example is the two- loop contribution to the anomalous magnetic moment.) They are pure short- 
distance corrections and arise from loop momenta of order nie and above. In the effective field theory 
approach they are included as a finite renormalization of the coefficients of the NRQED operators. 
These non-recoil corrections were first evaluated numerically a certain time ago [ p5| . More recently, 
they were calculated analytically by two independent groups ^] who agreed with each other but 
disagreed slightly with the numerical result (by about 0.10 MHz). The number we quote in the table 
is based on the analytical expression. The error of the result at order ruea^ in row 1 is of the level of 
a few 0.01 MHz and not indicated explicitly. It comes from the uncertainties in the Rydberg constant 
(Ref. (3^]) and in a (Ref. |^2|). The errors given in rows 5 and 7 are of numerical origin. For all 
other contributions the errors are negligible. The uncertainties due to the ignorance of the remaining 
meCK^ In and mgO^ contributions are not taken into account in the summed results. 

Except for the recoil corrections, all the quoted results are by now well established. There is 
still some controversy concerning the recoil corrections for which three different results can be found 
in the literature. The first calculation was performed by Caswell and Lepage in their seminal paper 
on NRQED ||T^. Recently, new calculations were performed by three different groups, using different 
techniques. First, Pachucki |14], using an effective field theory approach in coordinate space and a 



different regularization scheme, obtained a result differing significantly from the one of Caswell and 
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Order 



Specification 



Analytical/ 
numerical 



Contr. in MHz Refs. 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



meor 
nieOC' 



non-recoil 

recoil (Caswell-Lepage) 

recoil (Pachucki, Czarnecki et at) 

recoil (Adkins-Sapirstein) 

1- photon annihilation 

2- photon annihilation 

3- photon annihilation 



a 
a 
a 
a 
n 
a 
n 
a 
a 
a 
a 



204386.7 
-1005.5 
19.1 
-10.43 
3.1(6) 
7.02 
1.32(7) 
-2.34 
-0.61 
-0.97 
-0.92 



m 



35, 30, 16l + [B| 



m 

this work, O 



§0 



Sum (Caswell-Lepage) 

Sum (Pachucki, Czarnecki et al 

Sum (Adkins-Sapirstein) 

Experiment 



203 388.1(6) 
203 392.1 
203 386.4 

203 389.1(7) g 
203 387.5(1.6) g, § 



Table 1: Summary of theoretical calculations to the hyperfine splitting and most the most 
recent experimental measurements. 



Lepage. Then the Bethe-Salpeter formalism was employed by Adkins and Sapirstein ||T^, yielding 
yet another result. Finally, Czarnecki, Melnikov and Yelkhovsky ||l^ performed the calculation in 
momentum space with dimensional regularization. They obtained an analytical expression which 
agrees with Pachucki's numerical result. The number quoted in the table (row 6) corresponds to their 
analytical expression. 

Comparing with the most recent experimental measurement from Ref. Q, and ignoring re- 
maining theoretical uncertainties, the result containing the Caswell-Lepage calculation for the recoil 
contributions leads to an agreement between theory and experiment (VFth — W^cx = —1.0(1.0) MHz), 
whereas the prediction based on the result by Pachucki and Czarnecki et al leads to a discrepancy of 
more than four standard deviations (VFth — W^x = 3.0(0.7) MHz). On the other hand, the Adkins- 
Sapirstein result differs by slightly less than four standard deviations but, in contradistinction with 
Pachucki's result, it lies below the measured value: (PFth — ^cx = —2.7(0.7) MHz). Recently, the 
NRQED calculation [p| has been repeated by one of us (P.L.), in collaboration with R. Hill. Pre- 
liminary results of this calculation are in contradiction with the original NRQED calculation and also 
in agreement with Pachucki and Czarnecki et al. Hopefully, the theoretical situation will soon get 
settled. If the result of Pachucki and Czarnecki et al is indeed confirmed, the significant discrepancy 
with experiment will have to be addressed. 

We note again that we have not included any estimates about the remaining theoretical 
uncertainties in the considerations given above. The next uncalculated corrections are of order 
rrii.a'^ lna~^ ~ 0.7 MHz and could significantly influence the comparison between theory and ex- 
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periment. However, the coefficient of tlie log corrections are usually much smaller than 1, and we 
therefore believe that a contribution of 1 MHz from the higher order corrections is probably a conser- 
vative estimate. In this case, the discrepancy between theory and experiment remains unexplained. 
Clearly, further work on positronium calculations is necessary. 



8 Summary 

We have provided the details of the NRQED calculation of the 0{m(,a^) contribution to the positro- 
nium ground state hyperfine splitting due to single photon annihilation reported in an earlier paper. 
The counting rules needed to this order have been explained in detail and a discussion on some of the 
issues related to the use of an explicit cutoff on the momentum integrals has been given. We have 
provided a list of integrals useful for the evaluation of non-relativistic scattering diagrams. Our result 
completes the 0{m(,a^) calculation of the ground state hyperfine splitting and permits a comparison 
between theory and experiment at the level of 1 MHz. A comparison with the most recent exper- 
imental measurement underlines the need for more theoretical work concerning the C(mea^) recoil 
corrections and higher order contributions. 
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A Useful Integrals for the Matching Calculation 

In this appendix we present a set of integrals that has been useful for the matching calculations carried 
out in Sec. |^. All integrals containing ultraviolet divergences are regulated by the cutoff A, where the 
relation A S> po is implied. Terms of order k > are discarded. As explained in Sec. § we 

have regulated all infrared divergences by using a small fictitious photon mass A. In the following we 
give the results for arbitrary values of A, and in an expansion for X <^ pQ discarding terms of order 
a'', k > 0. For the matching calculations presented in this work only the fully expanded results are 
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